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SUMMARY 


Fixed-base and two-body equations of motion for an Annular Momentum 
Control Device (AMCD) are presented. An AMCD consists of a spinning annular 
rim which is suspended by noncontacting magnetic bearings and powered b/ a 
noncontacting linear electromagnetic motor. The fixed-base equations are 
for a rigid AMCD rim suspended by magnetic bearings attached to a rigid 
fixed-base and the two-body equations are for a rigid AMCD rim suspended by 
magnetic bearings attached tc a rigid body spacecraft. The fixed-base 
equations are applicable to any potential ground based AMCD application such 
as energy storage. 


INTRODUCTION 


The Annular Momentum Control Device (AMCD) represents a new concept in 
momentum storage devices (ref. 1), and is proposed for the control and 
stabilization of spacecraft and other applications requiring momentum storage 
(ref. 2). Aside from the advantages of the AMCD as a momentum storage device, 
it appears that the concept may also have significant advantages as an energy 
storage device for Earth-based applications (refs. 3-5). The basic AMCD 
concept (shown schematically in fig. 1) consists of a spinning annular rim 
which is suspended by a minimum of three noncontacting magnetic bearing 
suspension stations and driven by a noncontact ing linear electromagnetic 
drive motor. The magnetic suspension stations consist of electromagnetic 
actuators which produce axial and radial control forces on the rim. The 
electromagnetic actuators are controlled by servos which use rim position 
sensors to obtain the required relative displacement information. The 
purpose of this paper is to present a set of fixed-base and two-body equations 
of motion for an AMCD which can be utilized for the analysis of suspension 
system control concepts for Earth-based .imi spacecraft applications. The 
fixed-base equations are for a rigid AMCD rim suspended by magnetic bearings 
attached to a rigid fixed-base and the two-bodv equations are for a rigid 
AMCD rim suspended by magnetic bearings .ittailu'd to a rigid body spacecraft. 


LIST OF SYMBOLS 


!'>} 

Initial displacement of the AMCD axis system from the 
inertially fixed 1^, 2^, 3^ axis system. 

I^sl 

External disturbance torque acting on the spacecraft 
written in the S system. 

I^a} 

Translational force on the AMCD rim written in the 
A system. 

KraI 

Radial force on the AMCD rim produced by a given 
magnetic bearing station. 

I^Rsl 

Radial force on the spacecraft produced by a given 
magnetic bearing station. 


Translational force on the spacecraft written in 
the S system. 

I^xaI 

Axial force on the AMCD rim produced by a given 
magnetic bearing station. 

l^xsl 

Axial force on the spacecraft produced by a given 
magnetic bearing station. 

^Rl’ ’^R2 

Radial forces on the AMCD rim resolved along the 
1^ and 2^ axes, respectively. 

Kt 

Total moment acting about the center of mass of the 
AMCD rim written in the A system. 

i=sf 

Reaction torque on the spacecraft caused by motion of 
the AMCD rim acting through the magnetic bearings 
written in the S system. 

1^1 

Relative displacement in inertial coordinates of the 
AMCD rim with respect to the spacecraft at a given 
bearing station. 

|®s| 

|g| transformed into spacecraft coordinates. 

\h\ 

Absolute angular momentum of the AMCD rim (AMCD 
components) 

l“sl 

Absolute angular momentum of the spacecraft 
(spacecraft components). 

)“aoI 

I^^A^ f^Aof 
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LIST OF SYMBOLS (cont'd) 



Moments of Inertia of iIk.> AMCD rim about principal 
body axes. 



Moments of Inertia of the spacecraft about principal 
body axes. 



Total moments acting on the spacecraft (spacecraft 
components) . 


I'aI 

I'dI 


Mass of the AMCD rim. 

Mass of the spacecraft. 

Vector, fixed in AMCD coordinates, which locates a 
point in the AMCD rim that coincides with the center 
point of a given bearing station. 

Displacement of the center of mass of the AMCD rim 
in inertial coordinates. 


r 

m 

I'aiI 

I'caI 

l''sil 


Length of |r^| (radius of the AMCD rim) . 

Spacecraft coordinates of a given bearing station. 

|r^[ in inertial coordinates. 

Position of the AMCD rim center of mass with respect 
to the origin of the inertial coordinate system. 

Position of the spacecraft center of mass with respect 
to the origin of the inertial coordinate system. 

|rg[ in inertial coordinates. 

Vector transformation matrix, AMCD componentc. to 
inertial components. 


[Tj,] Vector transformation matrix, spacecraft components 

to inertial components. 

Torques on the .\MCl) rJm produced by a given bearing 
station written in the A system. 

Torques on the spaeeer.il t produced a given bearing 
station written in t lu* system. 

Transformation matrix, AMCD hmly rates to AMCD Euler 
rates. 
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LIST OF SYMBOLS (cont'd) 


K\ 

I«aI 

l^st 

h\ 


|SaI 

I«aI 




Notation: 

[ ] 

[ ]" 

[]T 

1 ! 

L J 


AMCD rim translational rates written In the A system. 
Spacecraft translational rates written In the S system. 
AMCD Euler rates. 

AMCD Euler angles with respect to inertial axes 
developed with a 3» 2, 1 rotation sequence. 

Spacecraft Euler rates. 

Spacecraft Euler angles with respect inertial axes 
developed with 3, 2, 1 rotation sequence. 

Spin angle of axes l^g, 2^^ referenced to axes 
1 , 2 , respectively. 

t\ A 

Absolute angular velocity of the AMCD rim expressed in 
the A system. 

Angular velocity of the 1^, 2^, 3^ axis system written 
in the A system with respect to the I system. 

Angular velocity of the 1^, 2^^, 3^^ axis system with 
respect to the 1^, 2^, 3^ axis system, written in the 
AB system. 

Angular velocity of the Ig, 2g, 3g axis system with 
respect to the I system written in the S system. 

Orthogonal axis system. 


Rectangular matrix 


Inverse of 
Transpose of 
Column vector 
Row vector 


[] 


[] 
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LIST OF SYMBOLS (concl'd) 


A single dot over a symbol denotes a first derivative with respect to 
time. Double dots denote a second derivative. 


Subscripts: 
1, 2, 3 
A 
AB 

a. b, c 
S 


Orthogonal components along 1, 2, 3 axes, respectively. 
AMCD coordinate system. 

AMCD body-fixed coordinate system. 

Magnetic bearing station a, b, or c. 

Inertially fixed coordinate system. 

Spacecraft coordinate system. 


TWO-BODY EQUATIONS OF MOTION 
AMCD - Spacecrafc Equations 

This section presents the AMCD-Spacecraf t equations of motion in general 
terms. The coordinate systems used for the .^MCD are shown in figure 2. A 

set of orthogonal 2^, and 3^ body axes are fixed in the AMCD rim with 

the 3^ axis being the rim spin axis. In order to simplify the equations, a 

set of nonspinning 1 , 2 , and 3. axes are used to define the motion of the 

n A A 

AMCD rim with respect to inertial space. The 3^ axis is coincident with the 

3 axis and the remaining two axes lie in the plane defined by the 1 and 
AB Ad 

-AB 

A set of orthogonal Ig, 2g, and 3g body fixed axes defines the motion 

of the spacecraft with respect to inertial space. The initial alinement of 
the axis systems is shown in figure 3. The spacecraft axis system is 
initially alined with an orthogonal 1^, 2^, and 3^ system fixed in inertial 

space. The AMCD axis system is displaced from the I axis system by a 
distance D. The 1^, 2^, and 3^ axes are parallel to the 1^, 2^, and 3^ axes, 

respectively. 

The assumptions used for deriving the equations of motion are: 

1) The AMCD rim and spacecraft are rigid bodies, 2) the AMCD rim and 
spacecraft have negligible products of inertia, 3) the spin rate of the 
•\MCD rim is constant with respect to the 1^, 2^, 3^ axis systi'iu, and 4) the 
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spacecraft is nonspinning. 

AMCD .- Referring to figure 2, the 3^^ axis is an axis of symmetry so that 


^ABl " ^A1 “ ^AB2 " ^A2 " ^A 


The absolute (inertially referenced) angular velocity, of the 
can be written as ^ ^ 

I^aI = 1«a! * {«Ao! (2) 


a) 

AMCD rim 


where 


h\ 


is the absolute angular velocity of the 1., 2 , 3, system and 

i\ i\ A 


|^Ao| angular velocity of the 1^, 2^^, 3^^ system with respect to 

the 1^, 2^, 3^ system (for a general discussion of rigid body dynamics 

see ref. 6). The absolute angular momentum, |h^|» of the AMCD rim then c 
be written as ^ * 


where 


{"4 - [ia]I^a} 


(3) 


[^4 - 


0 0 


'a 0 


0 0 1 


A3j 


(4) 


|f.^o| given by (see figure 2) 

l^Ao}' * L“ 

From equations 2 and 3, can be written 


as 


1«a! ■ WI'^a! * k]l«Ao! 


Making the substitution 


results in 


1"a«! -['a]M 
1"a! ■ [v] K! ^ !“ao! 


(5) 


( 6 ) 


(7) 


( 8 ) 
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(9) 


I^Aoj given by (from eqn. 4 and 5) 

|«Aor-L® ° 'ai^saJ 

and is constant since is assumed consLnnt. The total moment 

acting about the center of mass of the AMCD rim and referenced to the AMCD 
axis system is then 


• Kl- 


where 




{=a} - !«a! " HI-aI 


( 10 ) 


is the skew symmetric cross product matrix defined by 


[«a] 


-Q 


A3 






A3 


-n 


~^A2 ^A1 


A2 

A1 

0 


(11) 


|g^| is the reaction torque on the i\MCD rim caused by motion of the 

spacecraft acting through the magnetic bearings. Differentiating equation 8 
results in 


■ [^a] 1-\! ^ l»Ao} 


( 12 ) 


I^Aoj zero since the spin rate of the AMCD is assumc-d constant. 

Substituting equations 8 and 12 into 10 and rearranging terms results in 

{«a! - [^Ar(!<=Al -[^-a] (['a]{«a[ M"ao1)) 


The AMCD rates are obtained by integrating equation 13. 
The AMCD Euler rates can be written 


1»a1 - [^a]1‘-a1 


(14) 


where 




*A2 

sin0, , 
A1 

tan0,. r.osO 

i 

COS0. , 
A1 

-“*""ai 

'a2 

sln0. , 
Al 

sec0,,, COS0 

i 


A1 


AlJ 


( 15 ) 


gives the transformation matrix from A.MCD body rates to AMCI) Euler rates for 
a 3, 2, 1 rotation sequence. 


The translational force on the AMCD rim can be written in AMCD 
coordinates as 


\\\- "A (W- [«a] {'a!) 


(16) 


where J F. I is the translational force of the magnetic bearings acting on the 

I AJ 

rim, m^ is the mass of the rim,jV^Jis the rim velocity in AMCD coordinates, 
is the cross product matrix defined above. Rearranging 16 results 


and 

in 


K] 


\\\ ■ 4 - ^a! - [«a] ^a! 


(17) 


The rim translational rates are obtained by integrating equation 17. The 


rim translational rates in inertial coordinates 


where 


egratlnj 

• 1^4’ 


become 


[M 



■ N (vaI 

(18) 

'a3^^®A2 

<"®A3"®A2®®Al-=®A3‘^®Al' 

"^®A3"®A2^®A1®-"®A3*®A1> 


'a3*^®A2 

<=®A3'®A2"®Al^"»A3‘^®Al' 

*®®A3®®A2'*'Ar‘^®A3®®Al* 

(19) 

1 

91 

CD 

> 

10 

^®A2*®A1 

"®A2"®A1 


vector 

transformation from AMCD 

coordinates to Inertial 



coordinates. Because of the size of the expression, sin has been shortened 
to s and cos shortened to c. The position of the rim center of mass with 
respect to the origin of the I coordinate system is, from figure 4, 


where 




is given by 


I'-caI - I"-.!! * i"! 

i^r ■ L“i “2 “4 


(20) 


( 21 ) 
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(23) 


(24) 


Spacecraft .- The momentum of the spacecraft In spacecraft coordinates is 

j«s) ■ 

The total moment, , about the center of mass of the spacecraft referred 
to spacecraft coordinates is 

1"4 ■ [*s] * W ['s] Ni 

where is the cross product matrix. |Mg| can be written 

|»sl ■ K} * 

where |Gg| are the reaction torques on the spacecraft caused by motion of 
the AMCD rim acting through the magnetic bearings and |Eg| are the external 
disturbance torques acting on the spacecraft, consists of torques due 

to the rotation of the AMCD rim plus torques due to the translational forces, 
produced by rim motion, which act through a point displaced from the center 
of mass of the spacecraft. These translational forces, cause a trans- 

lation of the spacecraft center of mass and can be expressed in spacecraft 
coordinates as 

^si- "s(|^s} h]j'’sl) 

where mg is the mass of the spacecraft and |Vg| is the spacecraft velocity. 
Rearranging 25 results in 


4 - [«s] H 


(26) 


The spacecraft translational rates are obtained by integrating equation 26. 
The spacecraft translational rates in inertial coordinates, 


where 


K] ^ 


j''csj ■ [^s] {''si 

w- 


become 

(27) 


f the same form as 


The position of the spacecraft 


center of .nas.s with respect to tlie I coordinate system i 


•’ l^csl 


(see fig. 4). 


Magnetic Be.nrlng Equations 

The magnetic bearing configuration considered is one that consists of 
three bearing stations spaced equidistant ly around the AMCD rim. The 
locations of tl>e hearings are shown in figure 5. Eacli bearing station, 
which is fixed to the spacecraft, consists of an axial and r.ullal bearing 
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segment and associated gap sensor for detecting the relative displacement of 
the AMCD rim. 


In obtaining the equations for the bearing forces and torques, the 
following aasumptlons are made: 1) The AMCD and spacecraft are free to move 

with respect to each other with the only restraining forces between them 
being the forces generated by the magnetic bearings, 2) the axial bearings 
produce no force components In the radial direction, and 3) the radial 
bearings produce no force components in the axial direction. 


Bearing Gap s.- The magnetic bearing forces are a function of the relttive 
displacement between the AMCD and spacecraft. In order to calculate this 
displacement at a given bearing station, define a vector, |rg|, which is 

fixed in spacecraft coordinates and which locates the center point of the 
bearing station under consideration. Next, define a vector, ♦ '"'hich 

is fixed in AMCD coordinates and which, when the coordinate systems are 
alined (see figure 5), locates a point in the rim that coincides with th 
center point of the bearing station under consideintlon. The vector jr^ 
in inertial coordinates becomes ' 


■ [^s] l^s} 


and 


h\ 


in inertial coordinates becomes 


^Al| “ W {"a1 


(28) 


(29) 


The displacement { g| » in inertial coordinates, of the i\MCD relative to the 
spacecraft bearing station under consideration then becomes 


"( Ki|) " ( {''cs} 


Since the magnetic bearing segments and associated sensors are fixed to the 
spacecraft, the equivalent change in bearing gaps at a given bearing station 
then can be approximated by transforming |g| into spacecraft coordinates and 
considering the axial and radial componencs at the center point of the station 
as illustrated in figure 6. In this figure the components of |g| for bearing 
"a" are shown in spacecraft coordinates. The axial gap change is 85^3 3 ^“^ 

the radif.l gap change is g^^^. gj^^ is obtained by resolving and gg^2 

along the radial line labeled r^. 

Bearing torques .- The torque on the spacecraft due to the forces .it a given 
bearing station will be 

\\s\ ■ h] t^xs! ^ [s] !^xs} 

where is the cross product matrix, Uj^gj axial (<•!••• on the 
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spacecraft produced by the bearing, and radial force on the 

spacecraft produced by the bearing. Sinn* Mie forces on the rim produced by 
the radial bearing segments are assumed in .irt In the 1^, 2 ^ plane, the 

torque on the AMCD rim produced by a glv«*n bearing station becomes 


1’''"*! ■ [^] {■'xa} 


where 


(32) 

is the axial force on the 


^r^J is the cross product matrix and 
rim produced by the bearing. The vectors |rg| and were defined 

previously. 


Linearized Equations 


In this section the equations of motion and magnetic bearing equations 
are expanded and linearized. In linearizing the equations, small angle 
assumptions are made and second order terms involving spacecraft and AMCD 
motion are neglected. Also, since the distance from the AMCD rim to a 
fixed point, i.e. , a magnetic bearing location, is unaffected by rotation 
of the AMCD axis system about the 3^ axis, motion about this axis is set 
to zero. 


AMCD rotational equation s. -Expanding the term 
products of rates results in 



and neglecting 



(32) 


Using equation 33, the expansion of equation 13 becomes 


( "ai 

1 

- "a2 »ao > 

J 

• 



^^A2 

^A3 

r' 

"a1 »A03> 

0 


(34) 


The AMCD body rates are obtained by integrating equation 34. 
angle and rate assumptions, tlie Euler rates given by equation 


( e 

j 


1 


A1 

^A2 



1 6 

1 1 


\ A3 


0 


By using small 
14 reduce to 


(35) 
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i 


AMCD translational aquations.- The term 

and will be neglected. Equation 17 then becomes 


consists of products of rates 


r-i 


^A1 

CM 

< 

•> 

■ 

^2 

1 • ' 





hi I 


Equation 18, again using small angle and rate assumptions, can be written as 


Finally, equation 20 becomes 


*’CA1 

*^dl 

+ 

‘‘i 

r 1 ■ ( 

CA2 

, 

*^d2 

+ 

■'2 

*^CA3 * 

*^d3 

+ 

■*3 


Spacecraft rotation .il e quations .- The spacecraft rotational equations can be 
expanded in the sat..e manner as the AMCD equations with equation 23 becoming 


• 

i ^ 

[ "SI 


j ^S2 


1 • 


' "S3 



S2I I ‘52(052 + Eg,) 

1 i.gj ) I 1 /Is3^^S3 ^S3^ 

The spacecraft Euler rates become 


r. 


Spacecraft translational equations .- The spa'-erruft translational equacirns 
can also be expanded In the same manner as the AMCD equations. Equation Zb 
becomes 


"si 


"si 

"S2 


^S2 

("S3 


^S3 


(41) 


and equation 27 becomes 


I ’^CSl 


1 • 

J ’^CS2 

■ "S2 

1 *^053 

'"S3 


(42) 


Bearing forces .- The vector |rg| for bearing station a (see fig. 5) can be 
written as ' ^ 


or 


j^Saf ■ [' 


(d, + r cos 60°) (d-+ r cos 30°) d 

i m z m 


3J 


l"s! 


- 1<‘*1 * <’*> “2 r„) djJ (.J1 


for bearing station b becomes 

T 


hb!" ' [<“1 * <"2 - 


|rg| for bearing station c becomes 

!"sc r ■ ii‘*i ■ '*2 ''3I 

|r^| for magnetic bearing stations a, b, and c, respectively, beomce 

i--Ab r ■ "J 


(45) 


The vector 


H 


' L- 


[) 


(46) 

(47) 

(48) 


spacecraft coord i n. it i-s h«ronu*s 


1 ‘^s} ■ [h] I 


(49) 
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From equdticns 30, 28, and 29, then 


|*s! ' I'caI ' [^s]^ [hf4) ^cs} * j's}) »«> 

using the same assumptions as used previously, T^aI 


'®A2"®S2^ 


(«S1-®A1> 


L®S2"®A2^ ^®Ar®Sl^ ^ 


Equation 50 for magnetic bearing "a" then becomes 


S*Sa! 


^'dl * ®S3^“^2 ' ®S2‘*3 


^'^d2‘‘'cS2* ■ ®S3*'^1 '*' 


<’^d3-'cS3> * «S2‘*r“si''2 * <®S2-®A2><‘*>'^m * 


(e^l-8s,) (V3/2)r; 


Similarly, for bearings b and c, respectively 


{*sb! 


<--dr-^CSl> «S3'<‘2-<VT/2)rJ - e„d 


m^ S2 3 


^’^d3“’^CS3^ ^S2‘^r®Sl'^2 ''' ^^S2"®A2^ " 

(e^l-9si)(yr/2)r^ 


I^Scj 


^'^dr'’csp '*' ®S3^3 ■ ®S2‘*3 
'’^d2-‘cS2> ■" ®s3<V‘‘l> * ^SlS 


<'d3-'cS3> " ®S2“l - ®Sl’*2 <®A2-®S2>"m 


14 


The radial gaps for bearings a, b, and c, respectively, become 

St 


Ra 

*Sal ^Sa2 

Rb 

*Sbl “ gg^2 

Sc 

" *Scl 


The torques on the spacecraft due to the axial and radial forces 
and c bearing stations, from equation 31, become 




(dj + (V5/2)r„)F^3^ - (V3/2)d3F33^ 
-('»)djF33^ + (-/3:'2)<il^RSa 


{\sb} " 


(d, - (>/372)rJFjj33 + (7T/2)d,F 


-(d, + (>i)r„)F^33 + (H)d3F33j, 


-<’“>‘‘2''RSb - (^3/2)djF33^ 


3 RSb 


jl j = 

( mSc 


^^2^X50 


"^‘^rV^^XSc " ‘*3^RSc 


‘^2^RSc 


Sinilarly, for torques on the AMCD, from equation 32, 


{^mAa| 


l^mAbl 






<^/'>VxAb 


(4)r F 
m XAb 


( 55 ) 


at the a, b. 


(56) 


(57) 


(58) 


(59) 


(60) 
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0 


l^mAc I 


The total torques on 
become 


■r F„. 
m XAc 


(61) 


the spacecraft and AMCD due to bearing forces then 


^s} ■ 

{\sa! 

{\sb} 

^ l^mSc! 

(62) 

- 

l^mAal 

{’^mAb! 

j^mAc) 

(63) 


FIXED-BASE EQUATIONS 


The fixed-base equations of motion are obtained from the AMCD-spacecraf t 
equations by making the spacecraft axis system a fixed reference system and 
setting the offset D to zero. The equations for the gaps become (equations 
52-54) . 




(^ 3 / 2 ) 


r b., 
m A1 


(.h) 


T 9 o 

m A2 



j’^d2 

(r^3 - (v/i-/2) - (b) 



dl 


d2 




(64) 


(65) 


( 66 ) 


The axial gaps, in terms of rim rotations and translation become 

(67) 


®Xal 


(i/3/2) -(h) r 1 

m m 


CD 

> 

00 

- 

-(»/572) r^ -(H) r 1 

m m 


®A2 

®Xc 


0 r 1 

m J 


|’^d3 
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The radial gaps in terms of translations become (from equation 55) 


*Ra 


1/2 

>/3/2 

*Rb 

- 

1/2 

-V3/2 

*Rc 

1 


-1 

0 


Cl 


( 68 ) 


The torques on the rim due to the axial bearings become (from eqn. 63) 


A1 


A2 

l^A3 


(7T/2)r„ -(y3/2)r^ 

m m 

-d/2)r^ -a/2)r 

m m 


m 

0 


XAa 


XAb 


XAc 


(69) 


and the radial forces on the rim resolved along the 1 and 2 axes become 
(see fig. 5) 


R1 


R2 


1/2 1/2 -1 

v^/2 v/3/2 0 


RAa 

'RAb 

‘rAc 


(70) 


From equations 3A, 35, and 38 the rim rotational and axial translation 
dynamics become 


(71) 


where 


"A1 


Id 


/h - 

=A1 - 

%3 

L ■ 


^A03 

'"'a ’^CA3 

^A 


F ■ F + 
A XAa 

^XAb^ "^XAC 


69 and 70 



[ ( 

-(✓3/2)r„, 

0 

B 

u 

CM 

1 


r 

m 

1 

1 

1 



(72) 


XAa 


XAb 


XAc 


(73) 
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Finally f from equations 36, 37, and 38 the radial translation dynamics 
become 



’^CAl 



"a 

'CA2 


^2 


where and Fj ^2 defined by equation 70, 


CONCLUDING REMARKS 


This paper has developed fixed base and two-body equations of motion for 
an Annular Momentum Control Device (AMCD). The fixed-base equations are for 
a rigid AMCD rim suspended by magnetic bearings attached to a rigid base and 
the two-body equations are for a rigid AMCD rim suspended by magnetic bearings 
attached to a rigid body spacecraft. The fixed-base equations are applicable 
to the analysis of any potential ground based AMCD application such as energy 
storage. 
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Figure 1.- Annul?*, r Momentum Control Device Concept. 
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Figure 3.- Initial AMCD-spacecraf t coordinate system alinement. 





Figure 3.- Lx>cation of AMCD magnetic bearings. 
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